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Microwave  Refractive  index  Structure  Function  Profiles 
(Cn?)  Measured  From  a  Small  Aircraft 


1.  INTRODUCTION 

In  September  1985  a  series  of  airplane  flights  were  performed  to  obtain  meas- 
urements  of  the  microwave  refractive  index  structure  function  (Cq)  in  and  around 
inversion  layers.  This  series  had  two  goals  associated  with  Cn2  measurement: 
one  was  to  provide  evidence  on  the  validity  of  the  measurements  themselves;  the 
other  was  to  provide  information  on  the  variation  of  Cn2  in  the  first  few  km  of 
the  atmosphere.  Of  primary  importance  was  the  thickness  and  magnitude  rf  any 
enhanced  layers.  The  basic  interest  in  the  measurements  of  Cn2  was  to  under¬ 
stand  the  performance  of  troposcatter  radio  performance.  The  magnitude  and 
shape  of  enhanced  layers  strongly  influence  the  performance  of  these  radios  and 
can  dictate  the  strategy  of  making  the  Cn2  measurements  themselves. 

The  tests  were  carried  out  at  Chatham.  Mass.  The  Cq2  measurements  were 
obtained  using  a  microwave  refractometer  with  the  sensing  cavity  mounted  to  the 
wing  strut  of  a  Cessna -172  aircraft.  Supporting  atmospheric  data  were  obtained 
from  a  portable  radiosonde  system  and  from  synoptic  radiosondes.  These  radio¬ 
sonde  flights  were  taken  at  the  National  Weather  Service  (NWS)  station  at 
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Chatham.  An  earlier  report*  compared  the  meaaurement  of  radio  ducta  by  the 
two  radiosonde  systems  and  the  airborne  ref rac  tome  ter. 


1  DESCRIPTION  OF  MEASUREMENT 
21  Mewuremefit  Tdakps 

A  time  sequence  of  refractive  index  variation  measurements  n(tj)  is  obtained 
with  a  microwave  refractometer  mounted  under  the  wing  of  an  aircraft.  This 
time  sequence  is  converted  to  a  power  spectral  density  function  4n  by  taking  the 
fast  Fourier  transform  (FFT)  and  dividing  by  the  frequency  separation 

*n(ft)  »  A(ft)1 2/Af  ,  (1) 

where  A(f)'s  are  the  amplitude  from  a  FFT  and  Af  is  the  width  between  frequen¬ 
cies  (fj). 

Assuming  the  Taylor  hypothesis  to  be  valid,  die  time  spectrum  is  then  con¬ 
verted  to  a  spatial  spectrum  by  dividing  by  the  average  velocity  of  die  air  relative 
to  the  bensor: 

*n(ki)  *  W/V  •  (2) * 

where  V  is  the  velocity  of  air  relative  to  the  sensor,  and  k  is  the  wavenumber. 
These  spectra  are  then  used  to  calculate  <1.2  according  to  the  -5/3  power 

2 

law  for  three-dimensional  turbulence  in  the  inertial  subrange  given  by  Tirtarski 

g 

and  Wyngaard  and  LeMone. 

Cn2  «  0.25  ^  (k)  k”5/3  .  (3) 

On  some  flights,  two  refractometers  were  used,  generating  two  time 
sequences  of  refractive  index  measurements  [  nj(tj).  n2^jH  •  The  coherence 
spectrum  and  phase  spectrum  for  tb°se  time  sequences  were  computed  using 


1.  Morrissey,  J. ,  Izumi,  Y. ,  and  Cote,  O.  (1986)  Intercomparisons  of  Radio¬ 

sondes  and  an  Airborne  Refractometer  for  Measuring  Radio  Ducts.  AFflL- 
TR-86-C143,  ADA175150. 

2.  Tartarski,  V.L.  (1961)  Wave  Propagation  in  a  Turbulent  Medium,  McGraw- 

Hill,  New  York. 

3.  Wyngaard,  J.  C. ,  and  LeMone,  M.  A.  (1980)  Behavior  of  the  refractive  index 

structure  parameter  in  the  entraining  convective  boundary  layer,  J,  Atmos. 

Sci.  37:1573-1585. 
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Coh(f) 


where  ^  ^  ii  the  coapectrum  of  end  n^, 

QD_  la  the  quadratic*  coapectrum  of  n1  and  n.,  and 
1  «  1  * 

4  n  i  _  are  the  power  spectra  of  n.  and  n_. 

1  1  1  * 


(4) 


2.2  FUght  Technique 

To  obtain  the  necessary  time  sequence  of  refractive  inch.*  measurements, 
die  sensing  cavity  of  a  microwave  refractometer  was  attached  to  the  wing  strut 
of  a  0172  aircraft  just  below  the  leading  edge  of  the  aircraft  wing  (Figure  1). 

The  aircraft  was  then  flown  at  various  altitudes,  maintaining  altitude,  air  speed, 
and  direction  as  constant  as  possible.  These  level  data  gathering  runs  were  at 
~  100  knots  true  air  speed  (TAS),  lasted  5-1/2  min  or  more,  and  covered  in 
excess  of  10  miles.  These  runs  were  from  about  5  miles  south  of  the  Chatham 
Weather  Service  Station  over  Monomoy  Island  to  about  5  miles  north  of  die 
station  beside  Orleans  (Figure  2). 

The  flight  levels  were  selected  for  several  reasons.  Usually,  one  level  was 
at  about  3  km  to  characterise  the  area  well  above  the  inversions,  if  any.  If  there 
were  inversions,  several  levels  were  selected  to  be  either  just  above,  in,  or 
just  below  the  inversion.  On  some  flights,  one  or  more  of  these  levels  were 
repeated.  Then,  when  time  permitted,  a  series  of  low  level  runs  were  per¬ 
formed  at  typically  15,  30,  60,  120,  and  240  m  altitude.  The  selection  of  die 
flight  levels  was  made  using  information  transmitted  to  the  aircraft  from  the 
radiosonde  runs.  In  addition,  each  flight  began  with  the  aircraft  spiralling  up 
from  the  surface  to  3  km,  and  data  from  this  ascent  was  also  used  in  selecting 
levels. 
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Figure  1.  Microwave  Cavity  Mounted  on  Strut  of  Cessna-172  Aircraft  as  Flown 
2.3  MoMurement  Sensora/Syatem 

H»e  primary  sensor  is  the  microwave  refractometer  previously  developed 

at  the  National  Telecommunications  and  Information  Sciences  (NTIS)  laboratory 
4 

in  Boulder,  Colo.  The  refractometer  is  a  single  cavity  Birnbaum  type.  Fig¬ 
ure  3  is  a  simplified  block  diagram  of  the  refractometer.  The  cavity  is  continu¬ 
ously  driven  by  a  signal  that  is  centered  on  the  resonant  frequency  of  the  cavity, 
has  a  10  kHz  sweep  frequency,  and  causes  a  deviation  of  ±40  kHz.  The  center 
frequency  is  made  to  track  the  resonant  frequency  of  the  cavity  by  feedback  from 
a  phase  sensitive  detector  that  offsets  the  voltage  controlled  oscillator  (VCXO). 
The  output  of  the  VCXO  is  averaged  and  scaled  to  produce  a  digital  output  that  is 
displayed  on  the  front  panel.  This  same  digital  output  is  passed  through  a 


4.  Thompson,  M.  C. ,  Marler,  F.E.,  and  Allen,  K.  C.  (1985)  Measurement  of 
the  microwave  structure  constant  profile,  IEEE  Trans.  Antennas  Propag. 
AF-28(No.  2):278-280.  - 
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Figure  2.  Map  Showing  Local  Area  of  Tests.  Chatham,  and  General  Flight  Path 


Figure  S.  Block  Diagram  of  Microwave  Refractometer 

digital 'to -analog  converter  and  represents  one  of  the  analog  outputs.  This 
channel  has  a  sensitivity  of  0. 01  V/n  units  and  an  averaging  time/update  rate  of 
1  sec  as  flown.  The  second  analog  channel  is  obtained  by  taking  the  analog  error 
voltage  and  passing  it  through  a  bandpass  filter  (l/lOO  to  10  cycles)  amplifier. 

This  channel,  as  flown,  had  a  maximum  sensitivity  of  1. 04  V/n  units  (1.28  V/n, 
refractometer  No.  2)  at  a  frequency  of  2  cycles. 

In  addition  to  the  refractometer,  the  air  temperature  was  measured  by  a 
bead  thermistor  mounted  just  above  the  refractometer  cavity,  and  pressure  was 
measured  inside  the  aircraft  by  a  Rosemount  Model  120 1F1  pressure  transducer. 

A  data  system  was  designed  for  the  measurements  program  (Figure  4). 

It  consists  of  a  Euromak  C  computer,  manufactured  by  Dr.  Weiss  GMBH,  that 
is  used  to  control  12  A/D  (12  bit)  channels.  Six  channels  are  programmed  to  be 
sampled  at  a  rate  of  two  times  per  second  and  six  channels  are  programmed  for 
64  times  per  second.  Only  two  of  the  fast  channels  were  used  in  this  series. 

The  system  is  under  the  control  of  an  operator  using  a  small  NEC  -  PC8201A 
computer.  The  operator  is  allowed  to  input  titles  and  observations  of  two  types; 
one  at  the  start  of  each  flight  or  new  tape,  and  another  in  between  each  data  run. 
The  first  is  used  to  give  the  location  and  general  observations  on  the  test  and  the 
atmosphere.  The  second  is  used  to  give  specific  data  and  observations  asso¬ 
ciated  with  the  upcoming  data  run  (i.  e. ,  air  speed,  altitude,  heading,  location, 
visibility,  and  bumpiness)  as  well  as  additions  and  comments  on  the  previous  run. 
The  data  were  stored  on  a  nine-track  magnetic  tape  recorder  (Kennedy  9832) 
during  'hese  tests,  but  the  system  has  since  been  modified  to  use  a  cassette  style 
tape  ret  order  (Kennedy  6455). 
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Figure  4.  Refractometers  With  Associated  Computers  and  Recorders  as 
Flown  in  Tests 


3.  DATA  PROCESSING 


Figure  5  symbolically  shows  the  data  processing  U3ed  to  obtain  each  esti¬ 
mate  of  Cn2.  Sixty-four  seconds  of  the  sampled  high  gain  channel  of  the  refrac- 
tometer  (4096  values)  are  detrended  and  converted  to  a  frequency  spectrum 
^n^V'  which  ^en  converted  to  a  power  density  spectrum.  This  step  required 
the  removal  of  the  frequency  dependency,  | ,  introduced  by  the  pre -whitening 
filter  in  the  refractometer.  This  frequency  spectrum  is  convex  ted  to  a  wave- 
number  spectrum,  <£n(IC),  by  dividing  by  the  true  air  speed  of  the  aircraft.  A 
sample  power  density  spectrum  from  the  2  5  September  flight  at  1618  EDT  at 
Chatham  is  shown  in  Figure  6.  This  spectrum  is  then  used  to  generate  the  esti- 
mate  of  Cn^  by  using  the  mean  value  of  Cn^  calculated  for  each  wavenumber 
corresponding  to  i  =  25  to  i  =  104.  These  represent  wavelengths  from  30  to 
130  m  at  an  air  upeed  of  51.4  m/sec  (100  knots),  which  is  representative  of  the 
speed  in  these  flights . 
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During  each  of  these  estimates,  84  sec  of  data,  the  aircraft  moves  about 
S  km.  Each  level  run  contained  five  such  segments.  Figure  7  shows  the 
results  from  the  flight  on  18  September  at  1528  EDT.  At  each  level  the  mean, 
usually  of  five  values,  tc  shown  and  the  standard  deviation  is  shown  as  bars.  The 
temperature  and  dew  point  profiles  from  the  radiosonde  closest  in  time  to  the 
aircraft  flight  are  also  shown. 

On  those  flights  where  two  refractometers  were  flown,  the  sampled  data  was 
processed  from  each  refractometer  separately  as  described  above.  This  pro¬ 
vides  time  series  of  dual  measurements  of  Cn2  where  the  refractometers  are 
4.25  m  apart.  Figure  8  shows  the  C^2  results  from  the  26  September  flight  at 
1434  EDT.  The  results  are  in  the  order  taken,  with  vertical  lines  used  to  sep¬ 
arate  areas  of  constant  altitude.  In  addition,  the  two  sets  of  sampled  analog  data 
were  used  to  generate  coherence  spectra  and  ph&Se spectra.  This  was  accom¬ 
plished  by  first  calculating  the  cospectra  and  quadrature  spectra  in  a  manner 
analogous  to  that  shown  in  Figure  5.  Figure  9  shows  a  typical  plot  of  these 
spectra  from  the  flight  of  25  September  at  1618  EDT. 


TEMPERATURE  CC) 


Figure  7.  Cn2  Profile  for  Flight  of  18  September  at  1525  EDT 
With  Temperature  (T)  and  Dew  Point  Temperature  (TD)  From 
Radiosonde  Flight  at  1510  EDT 
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1434-1534  EDT  26  SEP  85  CHATHAM  MR  f!2 


Figure  8.  Sequence  of  Dual  Cn2  Measurements  From  Flight  of 
26  September  at  1434  EDT 


Figure  9.  Sample  of  Coherence  and  Phase  Spectra  From 
Flight  of  26  September  at  1434  EDT 


4.  DISCUSSION 

Two  types  of  information  were  gained  from  the  program;  information  on  the 
measurement  technique;  and  data  on  Die  distribution  of  C^2  in,  the  lowest  couple 
of  km  of  Die  atmosphere. 

The  use  of  two  refrac  tome  ter  a  demonstrated  Die  high  degree  of  repeatability 
of  Die  technique.  This  can  be  seen  in  Figure  10*  which  shows  Die  digitised  volt¬ 
age,  sequences  from  a  couple  of  data  blocks  (512  samples.  8  sec)  from  Die  flight 
on  25  September.  This  high  degree  of  repeatability  is  also  evident  in  the  five 
figures  in  Appendix  A  Diat  show  the  mean  Cn2  values  from  two  refractometers 
(Figures  A7,  A9*  All.  A12,  and  A1S>.  These  figures  show  a  consistent  offset 
in  Die  two  measurements.  This  offset  is  not  unexpected  in  that  time  and  circum¬ 
stances  prevented  an  absolute  calibration  of  the  second  refractometer.  This 
agreement  is  also  evident  when  we  look  at  Die  individual  values  obtained 
during  each  level  portion  of  flight  (Figure  8).  Again,  the  offset  between  Die  re¬ 
fractometers  is  evident. 

The  final  way  that  Die  dual  refractometers  were  compared  was  by  coherence 
and  phase  spectra  of  the  two  signals.  Figure  9  shows  Die  signals  to  be  well 
correlated  at  wavelengths  longer  than  25  m  and  to  have  a  zero  degree  phase  shift 
between  refractometers.  At  shorter  wavelengDts,  the  correlation  deteriorates 
and  the  phase  shift  is  no  longer  zero  until,  at  wavelengths  of  less  Dian  5  m,  the 
picture  is  chaotic.  This  behavior  is  expected  for  two  sensors  separated  by 
4.  25  m.  Th«d  signals  should  be  well  correlated  with  zero  phase  shift  until  the 
wavelength  is  about  four  times  the  separation,  or  20  m. 

The  measured  values  of  Cn2  range  from  a  high  of  3  X  10  13  m  '  to  a  low  of 
5  X  10“ 17  m“2/3.  Of  more  importance  to  understanding  the  performance  of 
troposcatter  radios  is  the  large  observed  vertical  gradients.  One  of  the  largest 
gradients  can  be  seen  at  about  1600  m  on  the  second  flight  of  18  September 
(Figure  A2).  Cn2  goes  from  3  X  10_1*  to  9  X  10  17  m  2/  3,  a  ratio  of  300  to  1, 
in  a  vertical  separation  of  150  m. 

On  several  flights,  we  see  a  correspondence  between  the  typical  inversion, 
rising  air  temperature  and  falling  dew  point  temperature,  and  a  sharp  decrease 
in  Cn2  with  height  (Figures  A2,  A3,  A4,  A6,  A10,  and  A14).  This  decrease,  at 
times,  exactly  corresponds  to  the  decr  ease  in  dew  point,  while,  in  other  cases, 
it  is  displaced  vertically  upward  a  couple  of  hundred  meters.  Often  there  is  also 


5.  Priestley,  J.T. ,  and  Hill,  R.J.  (1985)  Measuring  high-frequency  humidity, 
temperature  and  radio  refractive  index  in  the  surface  layer,  J.  Atmos, 
and  Oceanic  Technol.  2:233-251. 
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Figure  10.  Sampled  Voltage  Sequences  From  Both  Refractometers 
From  Flight  of  25  September  at  1125  EDT 
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a  f harp  increase  tn  C^8  with  increasing  height  Just  below  the  inversion.  This 
results  in  an  enhanced  lu  *r  of  Ca8  right  at.  and  sometimes  slightly  above,  the 
inversion  where  die  ratio  inside  die  layer  to  outside  is  typically  from  40  to  100 
timid  (Figures  AS*  A8,  Aid*  and  A14). 

On  those  days  where  there  were  low  level  runs  below  200  m,  re  see  die 
enhanced  levels  of  Ch8  associated  with  surface  friction.  1  X  10“lS  m-2^3,  which 
usually  falls  off  to  about  1  X  10~*3  by  S00  m.  This  decrease  is  often  not 

monotonic,  but  shows  a  pronounced  increase  with  height  in  die  first  100  m.  This 
increase  with  height  may  be  an  artifact  of  data  processing,  since  the  altitude  is 
less  than  some  of  the  wavelengths  used  in  calculating  Cn8. 


5.  CONCLUSIONS 

The  results  provide  strong  evidence  that  these  measurements,  obtained  from 
a  refractometer  mounted  on  a  small  aircraft,  can  provide  useful  information  on 
die  small  scale  fluctuations  of  n.  This  can  be  used  in  the  study  of  the  refractive 
index  structure  function  that  influences  many  modern  radar  and  propagation 
systems. 

Several  cases  showed  layers  of  very  limited  vertical  extent  (<  300  m),  where 
Cq2  increased  by  40  to  300  times.  The  sharp  layers  were  usually  associated 
with  a  temperature  inversion  and  an  associated  sharp  decrease  in  water  vapor. 
The  shallowness  of  these  layers  might  cause  them  to  be  misrepresented  in  any 
experiment  with  a  vertical  sampling  interval  greater  than  200  m.  The  shape  and 
position  of  die  layers  could  be  critical  to  understanding  the  performance  of  var¬ 
ious  radars  (airborne  surveillance  radar b  and  synthetic  aperture  radars)  and 
communication  systems,  i.e. ,  the  shifting  of  the  layers  by  about  400  m  could 
increase/decrease  delay  spread  on  a  troposcatter  radio  link  by  several  times. 
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Figure  Al.  Measured  Profile  (18  Sept  0730)  With  T  and 
Tq  From  Radiosonde  Closest  in  Time 
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Figure  A2.  Measured  C,^  Profile  (18  Sept  1525)  With  T  and 
Tp  From  Radiosonde  Closest  in  Time 
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igure  AS.  Measured  Cn*  Profile  (19  Sept  0719)  With  T  and 
t)  From  Radiosonde  Closest  in  Time 
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Figure  A4.  Measured  CjP  Profile  (20  Sept  0741)  With  T  and 
Td  From  Radiosonde  Closest  in  Time 
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Figure  AS.  Measured  Profile  (20  Sept  1113)  With  T  and 
Tn  From  Radiosonde  Closest  in  Time 
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Figure  A6.  Measured  C Profile  (21  Sept  0747)  With  T  and 
To  From  Radiosonde  Closest  in  Time 
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Figure  A8.  Measured  Cn2  Profile  (25  Sept  0653)  With  T  and 
Td  From  Radiosonde  Closest  in  Time 
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Figure  A10.  Measured  Cn2  Profile  (25  Sept  1626)  With  T  and 
Tj)  From  Radiosonde  Closest  in  Time 
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